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Within the framework of lattice QCD a high statistics computation of the nucleon axial and tensor 
charges is given. Particular attention is paid to the chiral and continuum extrapolations. 



1. Introduction 

The axial and tensor charges of the nucleon 
are defined by {s^ — — to|,): 

= 2s^Aq{^i) (1) 

■niN 

with qi^lbq, T^"" = iqat'^-f^q. The chi- 

rality and charge conjugation of these opera- 
tors are (+, +) and (— , — ) respectively. The 
charges have a parton model interpretation: 

Aq{p) = / dx[{q^{x,fi) - qi{x,fi)}+ 
Jo 

{q^{x,fi) -q^{x,^i)}] 

Sq{n) = / dx[{q^{x,fi) - qTix,fi)}~ 
Jo 

{qj_{x,^i) -q-f-{x,n)}] (2) 

where x is the fraction of nucleon momentum 
carried by parton quark density g, (a;,/i) at 
scale fi, in scheme S. Aq is related to the low- 
est moment of the gi structure function and 
can be measured in (polarised) DIS, while hi 

*Talk given by R. Horsley at DIS99, DESY-Zeuthen. 



(for Sq) having a — chirality must be found in 
a reaction allowing the quarks to have a dif- 
ferent chirality, such as (polarised) Drell-Yan. 

Known for the charges is that Au — Ad — 
gA ~ 1.26, (from neutron decay) and in the 
heavy quark limit Au = 5u = 4/3, Ad = 6d = 
— 1/3. Indeed in the non-relativistic limit as 
fermions (ie here quarks) are eigenstates of 70 
then we expect Aq — Sq. 

2. The Lattice Approach 

The hypercubic discretisation of the Eu- 
clidean path integral for QCD with lattice 
spacing a and resulting Monte Carlo evalua- 
tion of the partition function allows, in princi- 
ple, a fundamental test of QCD. That being 
said, the lattice programme is rather like an 
experiment: careful account must be taken of 
error estimations and extrapolations. There 
are three limits to consider: 

1. The box size must be large enough so 
that finite size effects are small. Cur- 
rently sizes of ~ 2fm seem large enough 
(the N diameter is ~ 0.8fm). 
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Fi gure 1. The x-extrapolation for qa- 

2. The chiral limit, when the quark mass 
approaches zero. It is difficult to calcu- 
late quark propagators at quark masses 
much below the strange quark mass. For 
each a value, three (or four) heavier 
quark masses are used and a linear ex- 
trapolation is made to the chiral limit, 
see Fig. 0. (Within our precision, there 
is no difference to extrapolating to the 
u/ d quark mass or to the chiral limit.) 

3. The continuum limit, ^ (the lead- 
ing order discretisation effects of the ac- 
tion can be so arranged to be O(a^), 
[ |l|). We have performed simulations 
at three a values corresponding approx- 
imately to (0.093, 0.068, 0.051) fm, 
[fl-i = (2.12, 2.90, 3.85)GeV] using scale 
ro = 0.5fm, [r^^ = 395MeV], [|. 

Our lattice calculation of matrix elements is 
standard, see eg [ ^ . We only note here that 
we are working in the 'quenched approxima- 
tion' when the fermion determinant in the 
partition function is set equal to one. The 
nucleon matrix elements consist of two dia- 
grams: a quark line connected piece and a 
quark line disconnected piece {qldis). Only 
the quark line connected piece is calculated. 
Thus, strictly speaking, from eq. (||), we can 
only compute 5u — Sd and Au — Ad. However, 
due to the additional negative sign in eq. (^ 
for Sq, we might also expect to be able to de- 
termine (5m, dd separately. (Both the above 
described approximations allow considerable 
savings in computer time.) 

3. Renormalisation 

In general operators (or matrix elements) 
must be renormalised 0^(/i) = Z^{iJL)Obare 
(in scheme 5, at scale /i) before they can be 



compared with experimental results. For the 
axial current ensuring PC AC on the lattice 
gives Z_4 and has been calculated in [ IJ] . The 
tensor current is more complicated. It is first 
convenient to define a r(7i-operator, which is 
both scheme and scale independent by 

= Az^i^l)Tii^,) (3) 



with 



(4) 



exp 



dr,o 



The factor [AZ^'^]-'^ is plotted in Fig. |. 
Given T^s^ then from Fig. || and eq. (^ we 
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Figure 2. [AZ|/'^] ^ using the one and twoJpop 
expansions of the beta and gamma functions, [ BI in 

units of A*^'^ for rij = flavours. The relation to the 
scale ro is given in [ ^. So, eg, at = 1, 2, 5 GeV, 
[AZ|7^]-i = 0.959^, 0.91l|, 0.870| respectively. 

can thus find 7^^'^(/i) at any desired scale \i. 
For Z^^ only the result for one loop perturba- 
tion theory is known. A non-pcrturbative re- 
sult would be desirable, at present we shall use 
'tadpole improved' perturbation theory, which 
re-expands the perturbation series, [ ^, using 
a physical coupling constant and removes the 
lattice 'tadpole' diagrams. Our explicit proce- 
dure using g^^^ is described in [ |j. To check 
that this gives a plausible result, we calcu- 
late the p-tensor decay constant, which also 
requires the same renormalisation constant. 
This decay constant is defined by 

(0|r'"'|p(p, A)) = ^ {elf - elpn f^{f^) (5) 
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This gives for /p'"^^V™p the results shown in 
Fig. ^. BaU and Braun, [ ^ , using the sum rule 
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Figure 3. The ^ O limit for fp^'^^/mp. 

approach find a value for fp{pL — IGeV) of 
160± lOMeV (n/ = 3). Using Fig. |, this con- 
verts to f^'^^yrup = 0.217(14). From Fig. | 
we see that there is good agreement between 
the lattice and sum rule methods. 

4. Results 

We now present our results. In Fig. |^ we 
show gA = Au — Ad. There appears to be 
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Figure 4. The ^ O limit for qa and {6u - Sdys'. 
Also shown is the experimental result for qa- 

an gradient (in distinction to the other re- 
sults). It seems difficult to reach the phe- 
nomenological result. Perhaps this is due to 
a quenching effect. Also in this figure we 



show {5u — Sdy^^ . It would seem that we 
have a milder 0{a^) gradient and also that 
{5u ~ ddy^^ > gA- This might indicate that 
a non-relativistic description of spin structure 
for the 'quenched' nucleon is reasonable. Note, 
however, that we are well away from the non- 
relativistic limit of 5/3. To conclude, in the 
Table, we give our Uf — continuum results: 



gA — Au — Ad 
Am 
Ad 


1.14(3) 

0.889(29) +Aqgid^s 
-0.236(27) +AqgM^s 




0.213(6) 


{Su - Sdy3' 
{Suys' 
{Sdys^ 


1.21(4) 

0.980(30) 

-0.234(17) 
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